Solvable Model Research Articles

Numerical simulation of the coupled dynamics model of the projectile engraving process

AbstractThe projectile engraving process directly influences the projectile motion in‐bore and impacts the firing accuracy, firing safety, and barrel life of the gun. For this reason, attention has been focused on this research topic. To address the limitations of the “instantaneous engraving” hypothesis adopted in the classical interior ballistic theory, the VUAMP user subroutine, one of ABAQUS's secondary development interfaces, is utilized in this paper to realize the modeling and numerical simulation of a coupled dynamics model of the projectile engraving process. In addition to facilitating engineering applications, a polynomial fitting formula of the engraving resistance obtained by simulation is proposed and then used as a supplement to establish a closed and solvable interior ballistic model considering the projectile engraving process. By comparing with test data, the simulation accuracy of the coupled dynamics model is verified. Simulation results reveal that the engraving process takes 3.8 ms, accounting for 26% of the whole launch process, which takes 14.6 ms, demonstrating that the process is not instantaneous. The results of this paper can serve as a reference for future studies on the coupled solution of the projectile engraving process and interior ballistics of guns or gun‐like equipment.

Solvable model of driven matter with pinning

Solvable model of driven matter with pinning

Randomness, Integrability and Universality

In spring 2022, the Galileo Galilei Institute for Theoretical Physics hosted a seven-week Workshop on ‘Randomness, Integrability and Universality’. The Workshop addressed a series of questions on exactly solvable models of statistical mechanics, having numerous ties and overlaps with various problems in random matrix models, probability, representation theory, and combinatorics. Much recent progress in these areas exploits the underlying notion of quantum integrability. Here we report on the scientific motivations and background of this activity and on its main outputs.

Footprints of loop extrusion in statistics of intra-chromosomal distances: An analytically solvable model.

Active loop extrusion-the process of formation of dynamically growing chromatin loops due to the motor activity of DNA-binding protein complexes-is a firmly established mechanism responsible for chromatin spatial organization at different stages of a cell cycle in eukaryotes and bacteria. The theoretical insight into the effect of loop extrusion on the experimentally measured statistics of chromatin conformation can be gained with an appropriately chosen polymer model. Here, we consider the simplest analytically solvable model of an interphase chromosome, which is treated as an ideal chain with disorder of sufficiently sparse random loops whose conformations are sampled from the equilibrium ensemble. This framework allows us to arrive at the closed-form analytical expression for the mean-squared distance between pairs of genomic loci, which is valid beyond the one-loop approximation in diagrammatic representation. In addition, we analyze the loop-induced deviation of chain conformations from the Gaussian statistics by calculating kurtosis of probability density of the pairwise separation vector. The presented results suggest the possible ways of estimating the characteristics of the loop extrusion process based on the experimental data on the scale-dependent statistics of intra-chromosomal pair-wise distances.

Variational Tensor Wave Functions for the Interacting Quantum Spin Hall Phase.

The quantum spin hall (QSH) phase, also known as the 2D topological insulator, is characterized by protected helical edge modes arising from time reversal symmetry. While initially proposed as band insulators, this phase can also manifest in strongly correlated systems where conventional band theory fails. To overcome the challenge of simulating this phase in realistic correlated models, we propose a novel framework utilizing fermionic tensor network states. Our approach involves constructing a tensor representation of the fixed-point wave function based on an exact solvable model, enabling us to derive a set of tensor equations governing the transformation rules of local tensors under symmetry operations. These tensor equations lead to the anomalous edge theory, which provides a comprehensive description of the QSH phase. By solving these tensor equations, we obtain variational ansatz for the QSH phase, which we subsequently verify its topological properties through numerical calculations. This method serves as an initial step toward employing tensor algorithms to simulate the QSH phase in strongly correlated systems, opening new avenues for investigating and understanding topological phenomena in complex materials.

Solvable Model of Quantum-Darwinism-Encoding Transitions.

We propose a solvable model of quantum Darwinism to encoding transitions-abrupt changes in how quantum information spreads in a many-body system under unitary dynamics. We consider a random Clifford circuit on an expanding tree, whose input qubit is entangled with a reference. The model has a quantum Darwinism phase, where one classical bit of information about the reference can be retrieved from an arbitrarily small fraction of the output qubits, and an encoding phase where such retrieval is impossible. The two phases are separated by a mixed phase and two continuous transitions. We compare the exact result to a two-replica calculation. The latter yields a similar "annealed" phase diagram, which applies also to a model with Haar random unitaries. We relate our approach to measurement-induced phase transitions (MIPTs), by solving a modified model where an environment eavesdrops on an encoding system. It has a sharp MIPT only with full access to the environment.

Energy optimal scheduling strategy considering V2G characteristics of electric vehicle

in recent years, large-scale electric vehicles connected to the power grid has a significant impact on demand-side resources. In addition, the vehicle-to-grid technology of electric vehicles introduces more uncertainties, which brings great challenges to the optimal scheduling of integrated energy systems. To solve this problem, this paper fully considers the uncertainty and fluctuation of wind power. Firstly, the objective function of joint configuration is constructed to maximize the utilization of wind power and minimize the planned generation cost of the system. Secondly, the linearized power flow equations are used to represent the relations among the state variables, and the second-order cone relaxation method is applied to deal with the branch power flow constraints, which is transformed into a solvable mixed-integer second-order cone programming model. Finally, the validity of the model is verified by an example. The system achieves obvious peak-shaving and valley-filling effect and improves the economic benefit of the system.

Emergence of a diverse array of phases in an exactly solvable model

Emergence of a diverse array of phases in an exactly solvable model

Summation formulas generated by Hilbert space eigenproblem

We demonstrate that certain classes of Schlömilch-like infinite series and series that include generalized hypergeometric functions can be calculated in closed form starting from a simple quantum model of a particle trapped inside an infinite potential well and using principles of quantum mechanics. We provide a general framework based on the Hilbert space eigenproblem that can be applied to different exactly solvable quantum models. Obtaining series from normalization conditions in well-defined quantum problems secures their convergence.

The impact of memory on learning sequence-to-sequence tasks

The recent success of neural networks in natural language processing has drawn renewed attention to learning sequence-to-sequence (seq2seq) tasks. While there exists a rich literature that studies classification and regression tasks using solvable models of neural networks, seq2seq tasks have not yet been studied from this perspective. Here, we propose a simple model for a seq2seq task that has the advantage of providing explicit control over the degree of memory, or non-Markovianity, in the sequences—the stochastic switching-Ornstein–Uhlenbeck (SSOU) model. We introduce a measure of non-Markovianity to quantify the amount of memory in the sequences. For a minimal auto-regressive (AR) learning model trained on this task, we identify two learning regimes corresponding to distinct phases in the stationary state of the SSOU process. These phases emerge from the interplay between two different time scales that govern the sequence statistics. Moreover, we observe that while increasing the integration window of the AR model always improves performance, albeit with diminishing returns, increasing the non-Markovianity of the input sequences can improve or degrade its performance. Finally, we perform experiments with recurrent and convolutional neural networks that show that our observations carry over to more complicated neural network architectures.

Partial Lyα\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document} thermalization in an analytic nonlinear diffusion model

During recombination, the cosmic background radiation is disturbed, in particular, by Lyman-alpha emissions from neutral hydrogen. It is proposed to account for the subsequent time-dependent partial thermalization of the Lyα\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document} energy content in an analytically solvable nonlinear diffusion model. The amplitude of the partially thermalized and redshifted Lyα\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document} line is found to be too low to be visible in the cosmic microwave spectrum, in accordance with previous numerical models and Planck observations.

Solvable models for 2+1D quantum critical points: Loop soups of 1+1D conformal field theories

We construct a class of solvable models for 2+1D quantum critical points by attaching 1+1D conformal field theories (CFTs) to fluctuating domain walls forming a “loop soup”. Specifically, our local Hamiltonian attaches gapless spin chains to the domain walls of a triangular lattice Ising antiferromagnet. The macroscopic degeneracy between antiferromagnetic configurations is split by the Casimir energy of each decorating CFT, which is usually negative and thus favors a short loop phase with a finite gap. However, we found a set of 1D CFT Hamiltonians for which the Casimir energy is effectively positive, making it favorable for domain walls to coalesce into a single “snake” which is macroscopically long and thus hosts a CFT with a vanishing gap. The snake configurations are geometrical objects also known as fully-packed self-avoiding walks or Hamiltonian walks which are described by an \mathrm{O}O(n=0) loop ensemble with a non-unitary 2+0D CFT description. Combining this description with the 1+1D decoration CFT, we obtain a 2+1D theory with unusual critical exponents and entanglement properties. Regarding the latter, we show that the log contributions from the decoration CFTs conspire with the spatial distribution of loops crossing the entanglement cut to generate a “non-local area law”. Our predictions are verified by Monte Carlo simulations.

A solvable model for graph state decoherence dynamics

We present an exactly solvable toy model for the continuous dissipative dynamics of permutation-invariant graph states of NN qubits. Such states are locally equivalent to an NN-qubit Greenberger-Horne-Zeilinger (GHZ) state, a fundamental resource in many quantum information processing setups. We focus on the time evolution of the state governed by a Lindblad master equation with the three standard single-qubit jump operators, the Hamiltonian part being set to zero. Deriving analytic expressions for the expectation values of observables expanded in the Pauli basis at all times, we analyze the nontrivial intermediate-time dynamics. Using a numerical solver based on matrix product operators, we simulate the time evolution for systems with up to 64 qubits and verify a numerically exact agreement with the analytical results. We find that the evolution of the operator space entanglement entropy of a bipartition of the system manifests a plateau whose duration increases logarithmically with the number of qubits, whereas all Pauli-operator products have expectation values decaying at most in constant time.

Dual symplectic classical circuits: An exactly solvable model of many-body chaos

We propose a general exact method of calculating dynamical correlation functions in dual symplectic brick-wall circuits in one dimension. These are deterministic classical many-body dynamical systems which can be interpreted in terms of symplectic dynamics in two orthogonal (time and space) directions. In close analogy with quantum dual-unitary circuits, we prove that two-point dynamical correlation functions are non-vanishing only along the edges of the light cones. The dynamical correlations are exactly computable in terms of a one-site Markov transfer operator, which is generally of infinite dimensionality. We test our theory in a specific family of dual-symplectic circuits, describing the dynamics of a classical Floquet spin chain. Remarkably, expressing these models in the form of a composition of rotations leads to a transfer operator with a block diagonal form in the basis of spherical harmonics. This allows us to obtain analytical predictions for simple local observables. We demonstrate the validity of our theory by comparison with Monte Carlo simulations, displaying excellent agreement with the latter for a choice of observables.

Nonadiabatic transitions during a passage near a critical point.

The passage through a critical point of a many-body quantum system leads to abundant nonadiabatic excitations. Here, we explore a regime, in which the critical point is not crossed although the system is passing slowly very close to it. We show that the leading exponent for the excitation probability can then be obtained by standard arguments of the Dykhne formula, but the exponential prefactor is no longer simple and behaves as a power law on the characteristic transition rate. We derive this prefactor for the nonlinear Landau-Zener model by adjusting Dykhne's approach. Then, we introduce an exactly solvable model of the transition near a critical point in the Stark ladder. We derive the number of excitations for it without approximations and find qualitatively similar results for the excitation scaling.

Strongly clustered random graphs via triadic closure: An exactly solvable model.

Triadic closure, the formation of a connection between two nodes in a network sharing a common neighbor, is considered a fundamental mechanism determining the clustered nature of many real-world topologies. In this work we define a static triadic closure (STC) model for clustered networks, whereby starting from an arbitrary fixed backbone network, each triad is closed independently with a given probability. Assuming a locally treelike backbone we derive exact expressions for the expected number of various small, loopy motifs (triangles, 4-loops, diamonds, and 4-cliques) as a function of moments of the backbone degree distribution. In this way we determine how transitivity and its suitably defined generalizations for higher-order motifs depend on the heterogeneity of the original network, revealing the existence of transitions due to the interplay between topologically inequivalent triads in the network. Furthermore, under reasonable assumptions for the moments of the backbone network, we establish approximate relationships between motif densities, which we test in a large dataset of real-world networks. We find a good agreement, indicating that STC is a realistic mechanism for the generation of clustered networks, while remaining simple enough to be amenable to analytical treatment.

Reflectionless pseudospin-1 Dirac systems via Darboux transformation and flat band solutions

This manuscript explores the Darboux transformation employed in the construction of exactly solvable models for pseudospin-one particles described by the Dirac-type equation. We focus on the settings where a flat band of zero energy is present in the spectrum of the initial system. Using the flat band state as one of the seed solutions substantially improves the applicability of the Darboux transformation, for it becomes necessary to ensure the Hermiticy of the new Hamiltonians. This is illustrated explicitly in four examples, where we show that the new Hamiltonians can describe quasi-particles in Lieb lattice with inhomogeneous hopping amplitudes.

Persistence time bound for subdiffusion based on multidimensional thermodynamic uncertainty relation: Application to an analytically solvable model.

The thermodynamic uncertainty relation (TUR) is an inequality showing the tradeoff relationship between the relative fluctuation of current observables and thermodynamic costs. It is one of the most important results of stochastic thermodynamics. There are various applications for TUR, one of which is the recent finding of thermodynamic constraints on the time window in which anomalous diffusion of Brownian particles can occur, including subdiffusion and superdiffusion, which are slower and faster than normal diffusion, respectively. These constraints are quite nontrivial because they are not generally derived from the asymptotic normal-diffusive behavior of the anomalous diffusion itself. In this study, we applied multidimensional TUR to the subdiffusion of Brownian particles obeying multivariate Langevin dynamics with a translationally invariant Hamiltonian in equilibrium. Multidimensional TUR is an improved TUR that includes information on another observable in addition to the one currently being considered. The use of an additional observable yields tighter bounds on the current fluctuation than those obtained using TUR. As an example, we demonstrated our theory using the one-dimensional Rouse model, which is known as a simple and analytically tractable model of polymer chains. We demonstrated that we improved the bounds for the persistence time of subdiffusion of the Rouse model, which became tighter as a more correlated observable with the current was used.

The role of initial system–environment correlations with a spin environment

Open quantum systems are a subject of immense interest as their understanding is crucial in the implementation of modern quantum technologies. In the study of their dynamics, the role of the initial system–environment correlations is commonly ignored. In this work, to gain insights into the role of these correlations, we solve an exactly solvable model of a single two-level system interacting with a spin environment, with the initial system state prepared by a suitable unitary operation. By solving the dynamics exactly for arbitrary system–environment coupling strength while taking into account the initial system–environment correlations, we show that the effect of the initial correlations is, in general, very significant and nontrivial. To further highlight the importance of the initial system–environment correlations, we also extend our study to investigate the dynamics of the entanglement between two two-level systems interacting with a common spin environment.

Coulomb and Higgs phases of G2-manifolds

Ricci flat manifolds of special holonomy are a rich framework as models of the extra dimensions in string/M-theory. At special points in vacuum moduli space, special kinds of singularities occur and demand a physical interpretation. In this paper we show that the topologically distinct G2-holonomy manifolds arising from desingularisations of codimension four orbifold singularities due to Joyce and Karigiannis correspond physically to Coulomb and Higgs phases of four dimensional gauge theories. The results suggest generalisations of the Joyce-Karigiannis construction to arbitrary ADE-singularities and higher order twists which we explore in detail in explicitly solvable local models. These models allow us to derive an isomorphism between moduli spaces of Ricci flat metrics on these non-compact G2-manifolds and flat ADE-connections on compact flat 3-manifolds which we establish explicitly for SU(n).